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Abstract 

The deconfinement transition in SU(Ac) Yang-Mills is investigated by Monte Carlo simulations of the 
gauge theory discretised on a spacetime lattice. We present new results for 4 < Ai'c < 8 (in particular, 
for Nc = 5 and Nc = 7), which are analysed together with previously published results. The increased 
amount of data, the improved statistics and simulations closer to the continuum limit provide us with 
better control over systematic errors. After performing the thermodynamic limit, numerical results 
for the ratio of the critical temperature over the square root of the string tension ^/(T obtained 
on lattices with temporal extensions Nt = 5, 6, 7, 8 are extrapolated to the continuum limit. The 
continuum results at fixed Nc are then extrapolated to A^'c = oo. We find that our data are accurately 
described by the formula T^j \fo = 0.5949(17) + 0.458(18)/A'^. Possible systematic errors affecting our 
calculations are also discussed. 

Keywords: SU(A^c) Yang-Mills Theories, Large Nc limit, Deconfinement Transition, Lattice Gauge 
Theories. 



1. Introduction 

In recent years, various lattice studies have 
been performed for SU(A'c) gauge theories in the 
't Hooft limit (see e.g. Refs. jlj E] for a review). 
As a result, on the one hand the old idea that the 
bulk of the physics is shared between SU(3) and 
the simpler SU(oo) |3HS] has been confirmed; on 
the other hand, solid bases have been provided for 
gauge-string duality studies aiming at describing 
QCD-like theories (see e.g. Ref. [6] for an early 
review of this field). 

On the lattice, one of the main areas of activity 
has been the finite temperature regime [7H21]. In 
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particular, it has been shown that the deconfine- 
ment temperature can be determined with very 
good accuracy [THOl HTJ UHl El]- This suggests to 
use the deconfinement temperature as the phys- 
ical scale in large Nc limit studies of observables 
at fixed lattice spacing [251 - I29] . which are often 
a useful intermediate step before performing the 
continuum extrapolation. The main motivation of 
this work is to complement existing results on the 
deconfinement phase transition by providing the 
value of the (pseudo-) critical coupling at various 
temporal extensions for the gauge groups SU(5) 
and SU(7), which have not been investigated at 
finite temperature before. Some of these results 
have already been used in our study of glueball 
masses at large Nc at the critical coupling for a 
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temporal extension Nt of six lattice spacings |28] . 
There, the inclusion of the Nc = 5 and Nc = 7 
data allowed us to increase the precision of the 
large extrapolation of the masses. The results 
of that study suggest that knowing the critical 
coupling at several values of Nt also for iVc = 5, 7 
can improve large Nc extrapolations in the con- 
tinuum limit. 

In addition, the calculations presented in this 
work provide us with an opportunity to revisit 
the extrapolation of the critical temperature to 
Nc = oo in the continuum limit. Besides adding 
the new results to existing lattice data, we inves- 
tigate possible systematic errors. In particular, 
the finite size studies of Refs. [SI [13] have been 
performed on Nt = 5 lattices and the results have 
been used to perform the extrapolation at other 
values of Nt with fixed spatial size Ng. Although 
this procedure is well justified in principle, there 
is the danger that, since the critical coupling at 
Nt = 5 for Nc> Q is close to the bulk phase tran- 
sition, the obtained value for the coefficient of the 
leading correction in 1/N^ in the thermodynamic 
extrapolation is significantly affected by finite lat- 
tice spacing artefacts. Since the determination of 
that coefficient is performed at a unique Nt, if the 
proximity to the bulk phase is a problem, the de- 
termination of the critical coupling in the thermo- 
dynamic limit at larger A^^^, and as a consequence 
the continuum limit of the critical temperature, 
will be affected by a systematic error. In the same 
spirit, we have performed calculations at iV^ = 7 
for 4 < Nc < 8, so that the continuum limit can 
be obtained by extrapolating the numerical data 
for the four values of A^^ = 5, 6, 7, 8. This allows 
us to estimate the influence of the Af = 5 point on 
the obtained numerical value at zero lattice spac- 
ing, and hence to check whether the continuum 
limit is affected by lattice artefacts related to the 
bulk phase transition or to the use of too coarse 
a lattice spacing in the extrapolation procedure. 
Finally, we have performed high statistics calcu- 
lations in SU(8), in most cases consisting of at 
least one million measurements, in order for the 
system at criticality to perform at least 8 round 
trips [tunnellings) between the conflned and the 
deconflned phase. The same criterion in terms of 



number of tunnellings has been used to determine 
the statistics for all the new simulations discussed 
in this paper. This should remove any bias related 
to a possible loss of ergodicity in the critical re- 
gion. 

The rest of this letter is organised as follows. 
In Sect.[2]we describe the system under investiga- 
tion, deflne the observables we study and provide 
numerical results for an estimator of the coupling 
at which the deconflnement phase transition takes 
place at flxed Nc , Nt and A^^. Sect. |3] deals with 
the thermodynamic limit of the critical couplings. 
Sect. |4]is devoted to the continuum extrapolation 
of the critical temperature in units of the string 
tension at flxed Nc- The large Nc limit of the 
latter quantity is discussed in Sect. [5j Finally, 
in Sect. |6] we summarise our results and briefly 
discuss possible future directions. 

2. The phase transition 

Our calculation follows the method exposed 
e.g. in Refs. [71 |8], which we will briefly sum- 
marise below. We consider a SU(A'c) gauge theory 
described by the Wilson action 

S = (3J2(^- JTc^'^ Tr {U,^{t))] , (1) 

where Ufj^yii) is the path-ordered product of the 
links U^{i) G SU(A'c) around the lattice plaquette 
identifled by the point i and the directions /i and 
p. /3 is deflned as (3 = 2Nc/gl, with the bare 
gauge coupling. The flnite temperature regime is 
realised by considering the system on a lattice of 
volume Ll X Lt, where Lg = aNg and Lt = aNt, a 
being the lattice spacing, with Nt Ng. Periodic 
boundary conditions are imposed in all directions. 
The temperature of the system is then given by 
T = 1/Lt and for fixed At it becomes a function 
of (3 only, through the dependence of the lattice 
spacing a on the gauge coupling. To find the value 
of the coupling at which the deconfining transition 
takes place, at fixed Nc, we compute the spatial 
average of the temporal Polyakov loop 



where x and t are the components of the Eu- 
chdean four-vector i respectively in the spatial 
and in the temporal directions (the latter corre- 
sponding to the dimension of size Lj). The decon- 
finement phase transition can be seen as a tran- 
sition from the phase symmetric under the center 
symmetry Z^r (the confined phase), to the phase 
in which this symmetry is spontaneously broken. 
Ip is the order parameter of the deconfinement 
phase transition. In addition, we study the four- 
volume average of the plaquette Up, defined as 

= E Tr {U,A^)) ■ (3) 

At fixed volume, we define the coupling /3c{Ns, Nt) 
corresponding to the deconfinement temperature 
by looking at the peak of the susceptibility of the 
modulus of Ip-. 

XI = N!m')-{\W) ■ (4) 

Our calculation is meant to complement the re- 
sults already present in the literature [3 El [131 EI] ■ 
Since calculations at larger than 8 become 
quite expensive |2T], we focused our attention to 
lower Nc- For Nc < 4, very precise results are al- 
ready available. For Nc = 4: and Nc = Q, previous 
studies already attained a good level of precision 
at Nt = 5, 6, 8. We hence studied the case Nt = 7, 
which has not been investigated before. The ad- 
dition of these calculations helps to improve the 
extrapolation of the corresponding /3c to the con- 
tinuum limit {Nt = oo). Calculations at Nc = 5,7 
have not been performed before. Hence, most of 
our numerical effort is devoted to the determina- 
tion of the critical temperature Tc for SU(5) and 
SU(7). Finally, for completeness, we have per- 
formed a new, high statistics numerical investiga- 
tion of SU(8) (comparable to that of [21]), which 
has enabled us to perform a more robust large Nc 
extrapolation. 

SU(A^c) gauge theories for Nc> 5 have a bulk 
phase transition at some value /3b of the coupling 
constant. The continuum physics is realised for 
/3 > 13b. For SU(5), /3b ~ 16.655 [13]. We 
have determined /3b for SU(7), which turns out 
to be around 33.246. The request that the sys- 
tem at criticality be in the continuum regime is 



A T 

Nc 


A T 

Nt 


A T 

Ns 


AT 1 r\ 
A^meas X 10^ 


4 


7 


22 


12 


5 


5 


8,10,12,14,16 


3 




6 


14 


2 




7 


16 


6 




8 


18 


5 


6 


7 


20 


10 


7 


5 


8,10,11,12 


2 




6 


10,12 


3 




7 


11 1 9 1 1 /I 







8 


14 


6 


8 


5 


7,8,10,11 


12 




6 


10 


6 




7 


12 


10 




8 


14 


10 



Table 1: Simulated volumes Ng for different gauge groups 
SU(A*'c) and different temporal lengths Nt. An approx- 
imate counting of the total number of measurements on 
each lattice is also shown. The results on these new lat- 
tices complement and improve the study of Ref. |13j . 

fulfilled a Nt > 5. As expected, this is the same 
bound on Nt already found for A^^ = 4, 6, 8. At 
fixed Ac, Ng and A^ we have computed xi (s^s 
Eq. (|4])) for about 10 /3s in the critical region, 
in a range that covers both the confined and the 
deconfined phase. For each calculation, we have 
used a combination of overrelaxation and heath- 
bath updates with ratio 4:1. The number of sweeps 
has been chosen in such a way that at least eight 
tunnellings were observed. In fact, in most of the 
cases we observed 12-15 tunnellings for the largest 
lattices. The statistics for each gauge group, spa- 
tial and temporal sizes is reported in Tab. [1} 

Always at fixed Ac, A's and Nt, using the points 
that are closer to the critical /3 (typically five or 
six), we have reweighted xi using the Ferrenberg- 
Swendsen procedure [30], as illustrated for one 
set of parameters in Fig. [TJ In particular, the 
reweighted data have been generated for several 
different bootstrap samples of the original simula- 
tions in order to give an unbiased estimate of the 
statistical error. One /3 value corresponding to the 
location of the maximum is chosen for each boot- 
strap sample (which is sufficiently well-behaved 
for a unique choice to be made). The central 
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Figure 1: Reweighted data of xi for SU(7) gauge theory 
on a 14"^ X 7 lattice. The filled circles are the measured 
data (with errors) and the unbroken line is the reweighted 
curve. The errors of the reweighted points (dashed lines) 
have been determined with a bootstrap (see details in the 
text). 



value and the standard deviation of the gaussian 
distribution of these (3 values determine our best 
estimate for the critical coupling f3c. The proce- 
dure described above allows us to avoid choosing 
the range for a quadratic fit approximating xi and 
therefore gives a more reliable and robust result. 

3. Thermodynamic limit 

For fixed and fixed Nf, the critical coupling 
/3c{Nt) is the thermodynamic limit of PdNs, Nf). 
Since for > 3 the phase transition is first or- 
der [71 IE], the extrapolation is performed accord- 
ing to the ansatz 



/3,(oo,iV,) = f3,{N,,Nt) + h{Nt) 
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where only the leading volume correction is taken 
into account and the value of the coefficient h{Nt) 
depends on the lattice spacing 
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The procedure to evaluate h{Nt) can be described 
as follows. At first we obtain the coefficient h{Nt = 
5) (corresponding to our largest lattice spacing) 
directly from a finite-size scaling (FSS) analysis 
using a wide range of volumes V = N^. We then 
use our own data for the beta function 
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ob- 
tained from the interpolation of the string tension 
over a large set of couplings in order to estimate 
h{Nt) at Nf = 6,7,8, with higher order correc- 
tions 0{a'^) accounted for by a 15% error increase 
on h{Nt) [13]. 

The determination of the beta function re- 
quires measuring zero-momentum correlators of 
Polyakov loops at zero temperature (torelons). In 
order to extract ami, the mass of the loop in lat- 
tice units, we look at the large separation expo- 
nential decay of those correlators, which is con- 
trolled by ami itself (see e.g. Ref. [25]). For 
SU(iVc) gauge groups with = 2,3,4,6,8, de- 
tailed measurements of torelon masses ami are 
already available on a wide range of coupling con- 
stants [25]. From the mass of such states ex- 
tracted using Polyakov loops of length Ng in units 
of the lattice spacing, we obtain the string tension 
ay/a by solving the equation 
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where the last two terms immediately derive from 
the effective theory describing the low-energy dy- 
namics of confining strings in the SU(A'^r;) the- 
ory [31]. If we keep a fixed to its continuum value, 
the numerical data for Oa/o" give us the variation 
of a as a function of /3. 

In Tab. [2] we summarise the string tension 
measured in high statistics simulations on large 
symmetric lattices L'^ at the reported values of 13. 
For the computation of a, we have used Eq. 
Since previous lattice calculations only used the 
leading correction —3^, we have verified that the 
insertion of the next-to-leading correction, whose 
universal character has been discovered only re- 
cently [211 E2], does not affect the numerical re- 
sults within errors. We have calculated for the 
first time the behaviour of the string tension in 
SU(5) and SU(7) function of the bare cou- 
pling. This allows us to have a precise estimate 
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Table 2: SU(5) and SU(7) string tensions on hypercubic 
lattices at the reported values of (3. The string tension 
is extracted from the mass of the lightest torelon state of 
length L. SU(4), SU(6) and SU(8) string tensions at the 
critical coupling for Nt = 7 are also shown. 

of the string tension at the couphngs correspond- 
ing to the deconfinement temperature for Nt = 
5,6,7,8. In addition to the above gauge groups, 
we obtained the string tension for SU(4), SU(6) 
and SU(8) in the neighbourhood of the critical 
couphng (3c for Nt = 7, needed for the improve- 
ment of the continuum extrapolation of the criti- 
cal temperature. 

For each value of N^ > 4, we interpolated the 
larger set of string tensions available to us using 
a polynomial function to obtain a{P)y/a. This in- 
terpolation includes a nested bootstrap sampling 
for a better error estimation and agrees very well 
with the fits of older data performed in Ref. [13] 
(for an alternative procedure, see also Ref. [5H]). 
Using a{l3)y/a and its derivative with respect to 
the coupling we obtain the infinite volume extrap- 
olation /3c(oo, Nt) at Nt = 6, 7, 8 shown in Tab. |3} 

As a final remark on the thermodynamic limit, 
we comment on the possible influence of the bulk 
phase transition. As it has been pointed out in 
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Table 3: SU(5), SU(7) and SU(8) values of the critical 
coupling for the corresponding temporal extent and in the 
infinite volume limit. Simulations are performed on x 
Nt lattices and the values f3c{oo,Nt) are obtained from 
Eq. ([5|. For SU(4) and SU(6) we also report our novel 
estimate of f3c at fixed Nt — 7. 

Ref. [2lj, at the lattice spacing corresponding to 
(3c for Nt = 5 the theory might still be affected by 
artefacts related to the nearby bulk phase transi- 
tion. In the SU(7) gauge theory at A^^ = 5, for 
which the bulk phase transition is very close to 
the finite temperature transition, we have explic- 
itly checked whether the determination of h[Nt) 
at that lattice spacing is significantly affected by 
lattice artefacts. To this purpose, we performed 
a FSS analysis on lattices with a larger tempo- 
ral extension {Nt = 7), where the deconfinement 
transition is pushed to a weaker couphng. Both 
procedures give infinite volume estimates for (3c 
that are compatible within the statistical uncer- 
tainty of our simulations. 

4. Continuum extrapolation 

Before we can take the continuum and the 
large Nc limit, the location of the deconfinement 
transition we found in the previous section needs 
to be translated into a physical temperature Tc. 
In addition, for the continuum limit, it proves con- 
venient to use dimensionless quantities, since at 
the leading order these have corrections that are 
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quadratic in the lattice spacing. Using the string 
tension a/o" to set the scale of pure gauge lattice 
simulations gives good control over systematic er- 
rors in the continuum extrapolation. For this rea- 
son, we study the continuum limit of the decon- 
finement temperature in units of the square root 
of CT, Tc/ yjo. This is determined for each iVc using 
only the leading O(a^) correction [8]: 



a=0 



+ be? a 



where at fixed lattice spacing 
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and 5 is a numerical coefficient of order one. We 
performed the continuum limit of T^j y/a accord- 
ma; to Eq. Q and using four different lattice spac- 
ings. The precision we achieve on the ratio T^/ -y/a 
is mostly determined by the precision of the string 
tension, since this latter quantity is affected by a 
relative error larger than that of /3c- 

The availability of an additional lattice spac- 
ing in the asymptotic scaling region for the con- 
tinuum extrapolation can help us identifying pos- 
sible systematic effects due to the inclusion of the 
coarsest point. In Fig. |2]we show the continuum 
limit of the deconfinement temperature for SU(5) 
and SU(7), which is the key original contribution 
of this work. Fits with and without the coars- 
est lattice point give compatible results for SU(7) 
with a psr d.o.f. above three disfavouring the 
former. For SU(5) the situation is similar, but the 
fit with all points has an acceptable x^- This is 
also true for N^. = 4, 6, 8. Our conclusion is that if 
there is any systematic effect in extrapolating the 
ratio Tc/ y/a to the a = hmit including points 
measured on A^^ = 5 lattices, this effect is signifi- 
cantly smaller than the statistical error. 

In Tab. H] we summarise the continuum limit 
values of Tc/ y/a that we used to obtain the SU(oo) 
limit. Only for A^^c = 7 we discard the coarsest 
lattice point in the continuum limit and show the 
result for the fit obtained using A^^ > 5. Since all 
the other gauge groups have well behaved extrap- 
olations with a low even including the Nt = 5 
point, in those cases we perform the fit using re- 
sults at all the available values of Nt. 




Figure 2: Continuum extrapolation of Tc/y/a for SU(5) 
and SU(7) according to Eq. (|9| and with different fitting 
ranges. Tlie continuum values are shown on the left. The 
dashed line, together with the open symbols, correspond 
to extrapolations without the coarsest lattice point. 



5. Large Nc extrapolation 

According to large A^c arguments, the large A'c 
limit of Tc/ a/o" can be expressed as a power series 
in 1/Nl 
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+ + oiNr' ,(11) 



where c is a numerical constant of order one. Our 
best fit for the large A^c deconfinement tempera- 
ture according to Eq. (11 ) is shown in Fig. [sj We 



extrapolate keeping only the 1/A^c correction to 
the planar limit, a procedure that has been shown 
to work very well down to A'c = 2 [13]. Includ- 
ing also SU(2) and SU(3) data from Ref. [H], we 
obtain 

TJ^/^ = 0.5949(17) + 0A58{18)/N^ , (12) 

with good x^/d.o.f. = 1.18. Discarding A"c = 2,3 
worsens the quality of the fit without changing the 
fitted parameters within the quoted error. Our 
value for the SU(oo) deconfinement temperature 
in units of the string tension is compatible with 
previous results reported in Ref. [13], but the rel- 
ative accuracy has increased approximately by a 
factor of 2. This is due to a better control over 
the continuum extrapolation for A"c > 4 and to 
the inclusion of SU(5) and SU(7) data. Note that 
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Table 4: Critical temperature in units of the string tension 
in the continuum limit for different gauge groups. The 
SU(2) and SU(3) values are taken from Ref. \TT . The large 
Nc extrapolation using all the reported values is shown in 
the bottom row. 



the more precise result is still compatible with the 
finite Nc value being accounted for by the leading 
correction only. 
In order to assess the robustness of our re- 
sult, we tested it against possible systematic er- 
rors in the continuum extrapolation. A different 
set of continuum values was created in the fol- 
lowing way: for each Nc > 4, results of contin- 
uum fits with and without the Nt = 5 point were 
merged together such that the error accounted for 
the whole possible range of values, while the mid- 
dle point of the error bar was taken as the central 
value. The estimates we obtained are fully com- 



patible with Eq. (12). Fitting the large Nc be- 
haviour of points obtained by extrapolating to the 
continuum limit results for Nf > 6 for all iV^ > 4 
(except for Nc = 5, where if we consider only 
points aX Nt > 5 the errors resulting from the 
continuum fit are anomalously small) also gives 
compatible results. 

6. Conclusions 

We have determined numerical values of the 
ratio Tc/ y/a for gauge groups SU(5) and SU(7). 
We have used the new data together with results 
for Nc = 2, 3, 4, 6, 8 already available in the liter- 
ature O El US] (supplemented with calculations 
at an additional lattice spacing for Nc = 4,6 and 
with calculations with increased statistics for Nc = 
8) to reanalyse the large Nc limit of this quantity. 
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Figure 3: Large Nc extrapolation of Td \fo using all data 
in Tab. [4] The dashed line corresponds to the fitted for- 
mula in Eq. ( 12 1. The inset is a close-up on the data. 



for which it turns out that only the leading l/A*"^ 
correction is needed to extrapolate finite A^^^ re- 
sults in the range 2 < Ai'c < 8. This had been 
already observed in previous simulations. We ob- 
tain an accurate large Nc limit that improves by 
a factor of two the precision of previous calcula- 
tions. At the same time, we investigated possible 
finite lattice spacing artefacts. Our analysis lead 
us to the conclusion that for 4 < Ac < 8, it is 
safe to extrapolate to the continuum limit from 
Nt = 5, as done in [71 [8]. 

In order to obtain a further noticeable im- 
provement, it is likely that gauge groups with 
Nc > 8 need to be investigated. However, since 
the strength of the first order deconfinement tran- 
sition grows with Nc, reliable Monte Carlo stud- 
ies of those systems will crucially require algo- 
rithms that mitigates substantially the exponen- 
tial suppression in the spatial volume of the tun- 
nelling rate between the confined and the decon- 
fined phases at criticality, like for instance the 
multicanonical algorithm [3l] . 
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